Optimal cloning of mixed Gaussian states 
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We construct the optimal 1 to 2 cloning transformation for the family of displaced thermal equi- 
librium states of a harmonic oscillator, with a fixed and known temperature. The transformation 
is Gaussian and it is optimal with respect to the figure of merit based on the joint output state 
and norm distance. The proof of the result is based on the equivalence between the optimal cloning 
problem and that of optimal amplification of Gaussian states which is then reduced to an optimiza- 
tion problem for diagonal states of a quantum oscillator. A key concept in finding the optimum is 
that of stochastic ordering which plays a similar role in the purely classical problem of Gaussian 
cloning. The result is then extended to the case of n to m cloning of mixed Gaussian states. 

PACS numbers: 



1. INTRODUCTION 

The no-cloning theorem states that quantum infor- 
mation cannot be copied, i.e. there exists no quantum 
device whose input is an arbitrarily prepared quantum 
system and the output consists of two quantum systems 
whose individual states coincide with that of the input 
0,0; Elm This and other quantum no-go theorems play 
an important role in quantum information theory and 
there exist deep connections with problems in quantum 
cryptography such as that of eavesdropping Q. With 
applications in mind, it is more interesting to derive a 
quantitative version of the theorem which says how good 
an approximate cloning machine can do, by providing 
lower bounds for the error made by any such machine. 

The quality of the approximate clones can be judged 
either locally, by comparing the state of each individ- 
ual clone with the input state, or globally by comparing 
the joint state of the approximate clones with that of 
independent perfect clones. Note that, because the no- 
cloning theorem requires that each individual system has 
the same marginal state as the input, it is the local qual- 
ity criterion which captures more of its flavor. However 
if we are interested in the joint state of the output then 
the global criterion is more useful as it takes into account 
the correlations between the systems. 

Before stating our cloning problem we would like to 
mention a few important results in this area and we refer 
to the review 5] for a more detailed discussion. The 
problem of universal cloning for finite dimension al p ure 
states was analyzed and solved in 0, 0, H, El EH 
IT3I Il4| . Interestingly, when the figure of merit is the 
supremum over all input states of the fidelity between 
the ideal and the approximate clones, it was shown that 
the same cloning machine is optimal from both the local 
and the global point of view 0, El • 

In the case of continuos variables systems the Gaus- 
sian states have received a special attention due to their 
importance in quantum optics, quantum communication 
and cryptography [l5|. Problems in quantum informa- 



tion such as entanglement measures ^ij and quantum 
channels 01 have- been partially solved by restricting to 
the framework of Gaussian states and operations. For 
coherent states th e op timal cloning problem has been in- 
vestigated in [l^, [H, under the restriction of Gaus- 
sian transformation. The question whether the optimal 
cloning map is indeed Gaussian has been answered pos- 
itively in the case of global figure of merit with fidelity, 
and negatively for the individual figure of merit [2l| . 

As noticed in 5], the area of optimal cloning for 
mixed states is virtually open partly due to the techni- 
cal difficulties compared with the pure state case. How- 
ever we should mention here the phenomenon of "super- 
broadcasting" H^,!^, [24| which allows not only perfect n 
to to (local) cloning of mixed states but even purification, 
if n is large enough and the input states are sufficiently 
mixed. This happens however at the expense of creating 
big correlations between the individual clones, just as in 
the case of classical copying. 

Quantum cloning shows some similarities to quantum 
state estimation, for example the pure state case is eas- 
ier than the mixed case in both contexts. Recently it 
has been shown that local cloning for pure states is 
asymptotically equivalent to estimation 25] . This pa- 
per makes another step in this direction by pointing out 
that global cloning has a natural statistical interpreta- 
tion. The statistics literature dedicated to the classical 
version of this problem pf3 | has been an inspiration for 
this paper and may prove to be useful in future quantum 
investigations. 

The problem which we investigate is that of optimal 1 
to 2 cloning of mixed Gaussian states using a global figure 
of merit. We show that the optimal doner is Gaussian 
and is similar to the optimal one for the pure state case. 
Our figure of merit is based on the norm distance rather 
than fidelity, the latter being more cumbersome to calcu- 
late in the case of mixed states. However the result holds 
as well with other figures of merit such as total variation 
distance between the distributions obtained by perform- 
ing quantum homodyne tomography measurements. 
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In quantum state estimation it has been shown |27| 
that the family of mixed Gaussian states appears as 
asymptotic limit of multiple mixed qubit states. Based 
on this result it can be proved [2^] that the problem of 
ri to 2ri global cloning of mixed qubit states is asymp- 
totically equivalent to that of 1 to 2 cloning of mixed 
Gaussian states which is addressed in this paper. 

In deriving our result we have transformed the optimal 
cloning problem into an optimal amplifying problem and 
then used covariance arguments to restrict the optimiza- 
tion to the set of mixed number states of the idler of a 
non-degenerate parametric amplifier p9l l3(| . The argu- 
ment leading to the conclusion that the optimal state of 
the idler is the vacuum, is based on the notion of stochas- 
tic ordering which is also used in deriving the solution to 
the classical problem of optimal Gaussian cloning. 

In Section|3|we extend the solution of the 1 to 2 cloning 
problem to the case of optimal n to m cloning of mixed 
Gaussian states. The transformation involves three steps: 
one first concentrates the n modes in one by means of a 
unitary Fourier transform, then amplifies this mode with 
a phase-insensitive linear amplifier with gain G = m/n, 
and finally the amplified state is distributed over the m 
output modes by using another Fourier transform with 
m — 1 ancillary modes prepared in a thermal equilibrium 
state identical to that of the input. 



2. CLONING OF MIXED GAUSSIAN STATES 

We consider the problem of optimal cloning for a family 
of Gaussian states of a quantum oscillator, namely the 
displaced thermal equilibrium states with a given, known 
temperature. Let a, a* be the creation and annihilation 
operators acting on the Hilbert space H. and satisfying 
the commutation relations [a, a*] = 1, and let 

oo 

$° := (1- S )5>»H, 

n=0 

be a thermal equilibrium state where < s < I is related 
to the temperature by s — e~ 13 , and \n) represent the 
Fock basis vectors of n-photons states. Let 

$ Q := D(a)<£°D{a)*, 

be the displaced thermal states where 

D(a) :— exp(aa* — eta), 

and consider the quantum statistical models: 

f:={* a :«eC}, e:={$"®$":«eC}. 

In the next Subsection we will give a statistical interpre- 
tation to the optimal figure of merit for cloning as a kind 
of gap (deficiency) between the less informative model J- 
and the more informative one Q . 



A. Figure of merit 

The aim of 1 to 2 global cloning is to transform the 
state $ a into <!>" ® $ Q without knowing a. This is how- 
ever impossible, and this fact has nothing to do with 
the quantum no-cloning theorem which is about local 
cloning. In fact the same phenomenon occurs in classi- 
cal statistics: given one Gaussian random variable whose 
distribution has unknown center, it is impossible to pro- 
duce two independent variables with the same distribu- 
tion. In both classical and quantum set-ups, if this was 
possible one could determine exactly the displacement by 
first cloning the state to an infinite number of indepen- 
dent states and then estimating the displacement using 
statistical methods. 

Thus we will try to perform an approximate cloning 
transformation which is optimal with respect to a given 
figure of merit. We consider a global criterion rather 
than a local, individual one. The classical version of this 
problem has been previously considered in mathematical 
statistics [2(J and we will adopt here the same termi- 
nology by defining the deficiency of the model T with 
respect to the model Q as 

Scion = 5(T,G) :=infsup||r($ a )-* a ®*°i|i 

where the infimum is taken over all possible cloning maps 
T : S(H) -> S(H) <g> S(H) with S(H) denoting the space 
of states (density matrices) on TL and T is a completely 
positive and trace preserving map. The norm one of a 
trace class operator is defined as ||r||i := Tr(|r|). We are 
looking for a map T = T opt satisfying 

F{T opt ) := sup ||T opt ($ Q ) - $ Q ® = Scion. 

This figure of merit is very natural from the statistical 
point of view |53 and can be related with the fidelity 
through the two sided inequalities (3l| 



-\\p-r\\l < 2 - 21^(^/2^1/2) < \\p-r\l,. 

Although the fidelity is a popular figure of merit, it is 
more difficult to handle in the case of mixed states. Note 
also that we do not take any average with respect to a 
prior distribution over the unknown parameter but just 
consider the doner which performs best with respect to 
all a, i.e. we are in a minimax framework as in |2l| . 

B. Cloning versus amplifying 

In the classical case the problem of optimal Gaussian 
cloning is equivalent to that of "amplifying" the location 
of the center of a Gaussian variable. We will show that 
this is also the case in the quantum setup by proving 
a fairly simple lemma allowing us to simplify the prob- 
lem and make the connection with the theory of linear 
amplifiers beautifully exposed in p9l |. 
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Let us start with the classical case and suppose that 
we draw a real number X from the normal distribution 
N(u, I) with unknown center ugl and fixed and known 
variance /. We would like to devise some statistical pro- 
cedure (for this purpose we may use an additional source 
of randomness) whose input is X and the output is a pair 
of independent clones (Y, Z) S ]R 2 , each having distribu- 
tion N(u,I). Let us assume for the moment that this 
can actually be done and note that by performing the 
invertible transformation 



(y,z)- (Qr + z)/V2, (y-z 



no statistical information is lost and moreover the two 
newly obtained terms are independent, (Y + Z)/y/2 has 
distribution N(>/2u, I) while (Y — Z)/y/2 does not con- 
tain any statistical information about u. Conversely, 
starting from an "amplified" version of X, that is a vari- 
able X with distribution N(\/2u,I), one can recover the 
independent clones by adding an subtracting an indepen- 
dent AT(0, 1) variable R: 



y= * + a Z= x-R 



V2 



V2 



In fact this is nothing else than the classical analogue 
of a 50-50 beamsplitter where Y, Z should be replaced 
by independent input modes carrying Gaussian states, 
and X, R can be seen as the output fields. The moral of 
this is that perfect cloning would be equivalent to perfect 
amplifying if any of them was possible, but in fact the two 
problems also are equivalent when we content ourselves 
with finding the optimal solution. We will prove this now 
in the quantum framework. Let A : S(Ti) — > S(Ti.) be a 
completely positive, trace preserving channel and define 
the figure of merit for y2 amplification 

F amp (A) :=8up||A($ a )-*^ , '|| 1 , 



and let S, 
merit. 



amp 



voIa F am p{A) be the optimal figure of 



Lemma 2.1 If T : S(H) m -> S(H) is an optimal 1 to 
2 cloning map then the map p i— ► Tri(£> o T(p)) is an 
optimal V2 amplifier, where B : S(H)® 2 -> S(H) m is 
the beamsplitter transformation which in the Heisenberg 
picture is given by the linear transformation 



( v/2 -y/2 



acting on the creation and annihilation operators of the 
two modes. Conversely, if A is an optimal amplifier, then 
the channel p i— > B~ 1 (A(p) ® $°) is optimal for cloning, 
and in particular 5 c i on — 5 amp . 

Proof. Let A be an optimal amplifier, i.e. F amp (A) = 



and T the corresponding cloning map, then 

F c lon{T) — 

= sup WB-^Aty ) ® *°) - B-\^ a (g) $°)||i 



sup \\A(<p a ) - = F amp (A) = S a 



Now let us suppose that there exists another clonig map 
T" with F c i n(T') < F c i on (T), then the corresponding 
amplifier A' satisfies 

F am JA') = sup ||Tn(i? o !"($<*)) - Tr^B^ ® ^ a ))\\ 1 



< \\T'(^ a ) - (g) $"||i = F clon {T') < F clon {T) = Samp- 

But this is in contradiction with the definition of the opti- 
mal figure of merit for amplification. A similar argument 
can be applied in the other direction. 

□ 



C. Covariance 

As in other statistical problems the search for an op- 
timal solution can be simplified if we can restrict the 
optimization set by means of a covariance argument. 

If the cloning map T has the property that if we first 
displace the input and then apply T, is equivalent to first 
applying T and then displacing the outputs by the same 
amount, then we say that T is (displacement) covariant: 



T{p) = D(-af z T(D(a)pD(-a))D(aY 



for all a £ C and p 6 S(H). By convexity of the 
distance we have 



F c i on (T) = sup||T($ a )-$ c 



> 



M ct ||T«($°) - $° ® $°||i > ||M Q T Q $° - $° <gi $°||i, 

where M a is the "mean with respect to a" , the analogue 
of averaging with respect to an invariant probability mea- 
sure for the case of compact spaces [2l|, |26| . Thus the 
mean M. a T a is at least as good as the initial channel T. 
There is a technical point here concerning the fact that 
M a r a may be singular as it is the case for example if T 
maps all states in a fixed one then M Q T Q = which is 
not trace preserving. A more detailed analysis 26] shows 
however that such cases can be excluded and one can re- 
strict attention to proper covariant and trace preserving 
channels. 

It can be shown [2l| that the general form of a covari- 
ant cloning map is given in the Heisenberg picture by the 
linear transformations between the input mode a,„ and 
the output modes ai and 02'. 



H = OA, + (K + hi) /y/2, 



a 2 = <2j. 



where b±, bi are two additional modes whose joint state 
determines the action of the cloning map. 
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The covariance property can be cast in the amplifier 
framework as well: a map A is a covariant amplifier if 

A(p) = D{-V2a)A(D{a)pD(-a))D{V2a), 

for all p € S(H) and a £ C. As shown in Lemma \2. II an 
optimal doner can be transformed into an optimal am- 
plifier by using a 50-50 beamsplitter to recombine the 2 
clones and then keeping one of the outgoing modes. For 
covariant cloning maps as described above in the rleisen- 
berg picture, this leads to the family of covariant ampli- 
fiers (see Appendix) with a amp := (ai + 0,2) 

a a mp = VZa in + b*, (2.1) 

where the mode 62 has been eliminated and the amplifier 
depends only on the state of the mode b := b\. Taking 
this into account, we will analyze the optimality problem 
in its formulation as optimal amplification. We will often 
use the fact that a particular covariant amplifier A is in 
one to one correspondence with a state r of the mode 
b as specified by the above linear transformation in the 
Heisenberg picture, and we emphasize this by writing 
A = A (r \ 

By using a further covariance argument we will show 
that the search for optimal amplifier can be restricted to 
states r which are mixtures of number states, i.e. states 
which are diagonal in the Fock vectors basis. Indeed for 
any displacement covariant amplifier A we have 

F amp (A) = sup \\A($ a ) - S^Hi = \\A(^ ) - 9%. 

a 

(2.2) 

Let Ub(0) = exp(i8Nb) be the phase transformation with 
Nb = b*b the number operator of the mode b, and de- 
fine similar phase transformations for the modes <Xj n and 
o-amp- The amplifier A is covariant with respect to phase 
transformations if 

A(p) = U am p(6)*A(U in {9)pU in (6)*)U amp (6) = A e (p). 

It is now easy to check that if A = A (r) then Ag = A^ 
where r g = Ub(—0)TUb(—0)*. Moreover, from (|2.2|l wc 
deduce that F amp (A) = F amp (Ag) because the state 3>° 
is invariant under phase transformations and thus 

F amp (A)= f d6\\A e ($ ) - 

> \\A($°)-<f>% = F arnp (A), 

where A = J AgdO = A^ and f is the phase averaged 
state r, i.e. a diagonal density matrix in the number 
operator eigenbasis. The rest of the paper deals with 
the problem of finding the optimal diagonal state for the 
mode b. 



D. Stochastic ordering 

Let us consider an arbitrary diagonal state r : = 
r «l n )( n l °f the mode b, and denote by p n = (1 — s)s n 



the coefficients of the thermal equilibrium state $° = 
^j n Pn\n){n\ of the a,i n mode. The state of the mode 
o-amp is itself diagonal and its coefficients can be written 
as Qi — 12k rPk T rtfk w hh fixed t\ k coefficients having a 
complicated combinatorial expression. The optimal am- 
plifier state T opt satisfies 

oc 

F amp (A^) = \\A (T °^(<P a ) - $ ||! = Y, \<h° pt ~Pi\ 

1=0 

00 

= inf F amp (A^) = inf V \qf - Pl \. 

T T Z * 

(=0 

The problem has been now reduced to the following "clas- 
sical" one: given a convex family T> := {q T : r £ S(H)} 
of discrete probability distributions on N and an addi- 
tional probability distribution p which does not belong 
to T>, find the closest point in D with respect to the || • ||i 
distance. In general such an optimization problem may 
not have an explicit solution but in our case the notion of 
stochastic ordering is a key tool in finding the optimum. 

Definition 2.2 Let p = {pi : I e N} and q = {q t : I S N} 

be two probability distributions over N. We say that p is 
stochastically smaller than q (p -< q) if 

m m 

/]pi - / Z®> Vm > 0. 

(=0 (=0 

The following Lemma is a key technical result which will 
allow us to identify the optimal amplifier map. 

Lemma 2.3 Assume that the mode a,i n is prepared in 
the thermal equilibrium state $° ; and the mode b in an 
arbitrary diagonal state t. Then the following stochastic 
ordering holds: q w ^< q T where q T is the distribution of 
the mode a amp defined in (|2.1|l and to = |0)(0| is the 
vacuum state. 

Proof. We will prove the result in two steps. First we 
show that the statement can be reduced to the case where 
the input state is the vacuum rather than a thermal equi- 
librium state <i> . Then we prove the lemma for the mode 
aj n in the vacuum state. 

In quantum optics the equation l|2.1|) describes a non- 
degenerate parametric amplifier [30] whose general input- 
output transformation has the form 

d(t) — cosh(x^)ci + sinh(xi)c2, 
C2(t) = cosh(xi)c 2 -I- sinh(x£)cj, 

where t represents the time and ^ is a susceptibil- 
ity constant. If both c\ and C2 modes are pre- 
pared in the vacuum state then each of the outputs 
separately will be in the thermal equilibrium state 
cosh(x£)~ 2 J2k tanh(xt) 2k \k)(k\. This means that we can 
consider that our input mode at n is one of the outputs of 
a parametric amplifier with tanh(x^) 2 = s. Thus 

din = cosh(x£)ci + sinh(x£)c2. 
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which together with (|2.1|l gives 

damp = v / 2cosh(xt)ci + V% sinh(xt) + b* 
= cosh(f)ci + sinh(i) (Tc* 2 + Rb*) 
— cosh(t)ci + sinh(t)c*, 

where cosh(i) = %/2cosh(x£), R = sinh(i) -1 and T — 
V 1 — R 2 . The right side of the last equation can be in- 
terpreted as follows: the modes C2 and b are combined 
using a beamsplitter with transmitivity T and one of the 
emerging beams denoted c is further used together with 
the mode ci, as inputs of a parametric amplifier char- 
acterized by the coefficient cosh(i). By hypothesis we 
assumed that the mode b is in state r, and by construc- 
tion the mode C2 is in the vacuum, thus the state of c is 
given by the well known binomial formula |32j 



for all to. With the notation 7 = T 2 we get 



1=0 

Jfe+m+1 \ ( fe ) 



1=0 



(1-7) 



fe+i 



kl 



1-7" 



1-7 



= 1 - 7 m+1 E( x - ^ k ~ r ( h+m+ v \ • 

r=0 



1_ 7 ™+I£ ( l_ 7 yy=- 



r=0 



Lemma 2.4 We have 



1 + S 



mo + 1 



□ 



. 2s mo+l 



The only property which we need here is that f is the 
vacuum state if and only if r is the vacuum state. In 
conclusion, by introducing the additional modes c\ and 
C2 we have transfered the "impurity" of the thermal equi- 
librium state from the mode aj„ to the mode c, and the 
stochastic ordering statement can be now reformulated 
in our original notations as follows: the mode ai n is pre- 
pared in the vacuum, and the mode b is prepared in a 
state f which is equal to the vacuum if and only if r is 
the vacuum. In addition, the relation H2.lt should be 
replaced by 

a-amp — cosh(i)<2i„ + sinhifo*. 

Under the assumption that aj„ is in the vacuum, we pro- 
ceed with the second step of the proof. Because stochas- 
tic ordering is preserved by taking convex combinations, 
we may assume without loss of generality that r = |fc)(fc| 
for k > 0. The following formula |30| | gives a computable 
expression of the output two-modes vector state of the 
amplifier 

if} = e r < n fc * e -9K^-+ fc * h + 1 ) e -ra ! „fc| 0j ^ 

1/2 



wheremo is the integer part o/log 2/(log(l + s) — log(2s)). 

Proof. Because both distributions are geometric, there 
exists an integer Too such that pi > qi for to < too and 
pi < qi for to > mo, and this proves the first equality. 
From the proof of Lemma l2~31 we can compute qf = (1 — 
s)s l where s = (1 + s)/2 and thus the integer too is given 
by the integer part of the log2/(log(l + s) — log(2s)). In 
conclusion 



lb-<?1i = 2 



1 + s 



mo+l 



2s 



-g{k- 



(=0 



l + k 



\l,l + k), 



where T = tanh(t) and e 9 — cosh(t). By tracing over the 
mode b we obtain the desired state of a amp 

1=0 i=a ^ ' 

. The relation q u ^ q T reduces to showing that 

tci m 

E%°>E^ 

;=o ;=o 



□ 

We arrive now to the main result of the paper. We 
will show that amplifier A T whose output is closest to the 
desired state state, is that corresponding to r = |0)(0|. 
Intuitively this happens because the "target" distribution 
q is geometrically decreasing and the closest to it in the 
family q T is the output which is the least "spread" . This 
intuition is cast into mathematics through the concept of 
stochastic ordering and the result of Lemma l2~3l 

Theorem 2.5 The state of the mode b for which the cor- 
responding amplifier map A is optimal is to = |0)(0|. In 
particular, the optimal amplifying and cloning maps are 
Gaussian. 

Proof. Define 



m a := max(m : ~^^Pi < a) 



1=0 



and T>(cl,t) = {D C N : T,ieD T i < °K for a11 a > °- 
Note that by Lemma 12.31 we have X/ISo 9 r — a a ^ 
r, and thus {0, 1, . . . , m a } S T> (a, r). Using the relation 
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\\p — q\\i = 2sup D ^2i eD (pi — qj) we obtain the chain of 
inequalities 

lb- <fl|i = 2 sup sup y2(pi-qj) 

a>0 DeV(a,r) 
m a m a 

> 2supJ2(Pl - if) ^ 2 SU P £ (Pi ~ if) 

1=0 1=0 
m 

= 2 S u P J2(pi-ii) = \\p-i u h- 



m>0 



1=0 



The first equality follows directly form the definition of 
2?(a, r). The following inequality restricts the suprcmum 
over all D £ T>(a, r) to one element {0, 1, ... , m a }. In the 
second inequality we replace the distribution q T by q u us- 
ing the stochastic ordering proved in Lemma 12.31 In the 
following equality we use the fact that both distributions 
p and q w are geometric (see also Lemma \l. 411 . 

As discussed in Section lTHl we can restrict to covariant 
amplifiers and the figure of merit in this case is simply 
||A T ($°) - $°||i = ||q T — thus the optimal ampli- 
fier is A" . Moreover, by the equivalence between optimal 
cloning and optimal amplification we also obtain the op- 
timal cloning map (see Lemma |2.1|) . 

□ 

Corollary 2.6 The optimal 1 to 2 cloning figure of merit 
is 

1 4- s\ mo + 1 
Salon = 2 [ -4- I - 2.S mo + 1 



with mo as in Lemma\ 
Proof. This follows from Theorem 12 . 51 and Lemma \l. 41 

□ 

E. Comparison with the classical case 

The derivation of our result on optimal quantum 
cloning is inspired by a similar one in the classical do- 
main |26j . In this subsection we comment on the optimal 
figures of merit in the two cases as function of the pa- 
rameter s = e~@ . 

It is well known that an arbitrary state p of a quantum 
harmonic oscillator has an alternative representation as 
a function W p : M 2 — ► R called the Wigner function. In 
the case of the family of displaced thermal equilibrium 
states the Wigner function is a two dimensional Gaussian 

M 



W go , Po (p,q) = -tanh(/3/2)e 



with variance 



-((9-9o) 2 + (p-po) 2 )tanh(/3/2) 



1 + 8 



and qo + ipo — V^a. We have shown that the best 
quantum amplifier produces a Gaussian state with s = 
(s + l)/2 or in terms of the variance 



3 + s 



2(1-*) 

which implies that for any s £ [0, 1) we have the relation 



V~ S = 2V S + 1 -, 



(2.3) 



2tanh(/3/2) 2(1 -s)' 



which indicates the least noisy amplification according to 
in the fundamental theorem for phase-insensitive ampli- 
fiers m. 

Let us consider now the classical problem of Gaussian 
cloning as discussed in the beginning of Subsection 12 Bl 
given a Gaussian random variable Igl 2 with distribu- 
tion N(u, I), we want to produce a pair (Y, Z) of indepen- 
dent clones of X . By using the equivalence between the 
cloning and the amplification problems, the task is equiv- 
alent to that of producing a variable X with distribution 
N(y/2u, I), and the optimal solution to this problem j2|| 
is simply to take X — \[2X\ We note that in the classi- 
cal case the variance of the output is always equal to the 
double of the variance of the input, while in the quan- 
tum case the output "noise" is always higher due to the 
unitarity conditions imposed by quantum mechanics |29| , 
and we recuperate the factor 2 in the high temperature 
limit (|2~3|) . 

In the classical case one can deduce by a simple scal- 
ing argument that the classical figure of merit does not 
depend on the variance of the Gaussian but only on the 
amplifying factor, and in our case it takes the value 1/2. 
As expected, the quantum figure of merit is larger than 
the classical one to which it converges in the limit of 
high temperature, s — ► 1. The upper line in Figure 
represents the optimal figure of merit 5 c i on — damp as 
function of s. An interesting feature of this function is 
that it appears to have discontinuities in the first deriva- 
tive precisely at the values of s for which the "crossing 
point " too makes a jump (see Lemma 12. 4|) . 

For comparison we have also plotted the norm one dis- 
tance between the corresponding Gaussian Wigner func- 
tions which does not seem to show any roughness. 



3. OPTIMAL n TO m CLONING OF MIXED 
GAUSSIAN STATES 

The results which we obtained for optimal 1 to 2 
cloning can be easily extended to the case of optimal 
cloning of n to to cloning of mixed Gaussian states. The 
idea is to first "concentrate" the state ($ Q )®" of the in- 
put into a single mode by means of a unitary transforma- 
tion followed by discarding the uninteresting n— 1 modes. 
Then, one amplifies the obtained state by a factor yjm/n 
(gain factor G — m/n) and distributes it using another 
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FIG. 1: Figures of merit for quantum and classical optimal 
cloning as function of s — e~' 3 . 



unitary transformation applied on the amplified mode to- 
gether with to — 1 additional ancillary modes prepared 
in state $°. The unitary transformations are the Fourier 
transforms 1241: 



a h . 



^ n— 1 

U n a k U* := —= y e 



In 



U m b k Ul 



1=0 



1=0 



at, 



(3.1) 



(3.2) 



where a k are the input modes and UmbkU^ are the out- 
put modes. The amplifying part is described by the co- 
variant map 



to — n , 



(3.3) 



where b is an additional mode prepared in a diagonal 
state r as in the 1 to 2 case, and 



^ n— 1 

:= U n a U* = -7= 22 



ai- 



v 1=0 

After amplification the second unitary transformation 
is performed on bo := a amp and the ancillary modes 
bi, . . . , b m -\ prepared in the state <&°. 

An obvious extension of Lemma \'2. II holds in this case 
as well, showing the equivalence of optimal cloning and 
optimal amplification. Similarly, Lemma |2. 31 and Theo- 
rem hold in general for any amplifying factor k > 1 
and we arrive to the conclusion that the optimal ampli- 
fier is given by the transformation described in equation 
H3.3(l with the idler mode b in the vacuum state. 

An interesting fact is that our transformations are sim- 
ilar to those of optimal n to to broadcasting [24| with the 
exception that in the last step of the procedure different 
ancillary states are used: the optimal state for broadcast- 
ing is the vacuum while for global cloning it is same ther- 
mal equilibrium state < I> which characterizes the family. 



We have constructed an optimal 1 to 2 cloning map 
for the family of displaced thermal equilibrium states of 
a fixed, known temperature. We have considered a global 
figure of merit based on the supremum over all displace- 
ments of the norm distance between the joint state of the 
approximate clones and that of the ideal ones. The op- 
timal doner is Gaussian and is similar with the optimal 
doner for coherent states with global figure of merit and 
consists of two operations. The amplification step uses 
a non-degenerate linear amplifier with idler prepared in 
the vacuum state. The cloning step uses a beamsplitter 
and another ancillary mode in thermal equilibrium state 
with the same temperature as the input. 

Computations which have not been included here in- 
dicate that the optimal cloning map remains unchanged 
under global figures of merit using different "distances" 
between states. 

The local version of the optimal cloning problem would 
probably lead to a non-Gaussian optimum as it is the case 
with coherent states |21| . 

The equivalence between cloning and amplifying can 
be extended to an arbitrary number n of input states 
and number of clones to, as well as the proof of the 
optimal amplifier. In the case n > 1, the first step is 
the concentration into one mode by means of a unitary 
Fourier transform, followed by amplification with gain 
factor G = m/n, and distribution into to output modes 
using another Fourier transform. 

Some other generalizations of the Gaussian cloning 
problem may be considered for future investigations, such 
as an arbitrary number of modes with larger families of 
Gaussian states. For example in the case of a family of 
thermal equilibrium states with unknown temperature, 
one may need to perform an additional estimation of the 
thermal states in the last step of the cloning which re- 
quires ancillary modes prepared in the equilibrium state 

Finally, the key ingredient in our proof was the notion 
of stochastic ordering which is worth investigating more 
closely in the context of quantum statistics. 
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APPENDIX: DISPLACEMENT COVARIANT 
AMPLIFIERS 

We give here a short proof of the fact that the dis- 
placement covariant amplifiers have the form Ij2.1|l . Let 
A : S(H) — > S(H) be a covariant amplifier such that 

A(p) = D(-V2a)A(D(a)pD(-a))D(V2a). 



8 



Then the dual A* : B(H) — > B(H) has a similar property, 
for all X e B(H), a EC 

A*(X) = D(a)A* (p(-V2a)XD(V2 a) ] D(— a). 

By choosing X = D((3) and using the Weyl relations 
D(a)D(/3) = exp(Hm(aj3))D(a +/?) we get A*(D(/3)) = 
c(@)D(y/2(3) for some scalar factor c(/3). Now, accord- 
ing to the Theorem 2.3 of 33] if A is trace preserving 



and completely positive the constant c(j3) is of the form 
c(a) = p{D(a)) where p is a state in S(TL). Thus 

A*(D(a)) = p{D(a))D(V2a). 

Now, it can be checked that if we start from (|2.1|) with the 
mode b prepared in state p then A describes the channel 
transformation from the input to the output mode. 
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